math-ph/0708.3687 
DAMTP-2007-78 
DCPT-07/47 



Multiplicity in Supersymmetric Spin Chains 

David Kagan"^ and Charles A. S. Young^ 

""DAMTP, Centre for Mathematical Sciences, University of Cambridge, 
Wilherforce Road, Cambridge CBS OWA, UK 

''Department of Mathematical Sciences, University of Durham 
South Road, Durham DHl 3LE, UK 
E-mail: d.kagan@danitp.cani.ac.uk, charles.young@durhani.ac.uk 

Abstract 

We discuss a simple procedure for obtaining new integrable spin chains from old by 
replacing each single state of the original model by some collection of states. This works 
whenever the Lax matrix of the chain has a certain form. The simplest example is the su(n) 
XX model. We apply the techniques of the nested algebraic Bethe ansatz to solve such 
systems, in the bosonic and supersymmetric cases. 

1 Introduction 

In [1], Maassarani and Mathieu introduced generalizations of the usual two-state XX spin chain, 
motivatated in part by the desire to construct versions of the Hubbard model with su(n) symmetry 
[21 [31 m [5] . Following [1] these su(n) XX models were further generalized in [6l[7] and it emerged 
that the underlying idea is rather simple. For a certain class of i?-matrices it is possible to 
introduce "multiplicity" . This means, roughly speaking, replacing each state of a given spin chain 
by a collection of states, all equivalent from the point of view of the spin chain interactions — 
we will make this more precise in section 2. In [7], this procedure was applied specifically to the 
i?- matrix of Am in the fundamental representation, but it is not limited to it. 

The main goal of the present work is to carry out a similar construction in the case of supersym- 
metric /^-matrices and spin chains. We first review, in section 2, the introduction of multiplicity 
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into XXZ-like integrable systems, focusing on identifying the features needed to make the con- 
struction work. Lax operators, monodromy, and transfer matrices are defined. In section 3 we 
examine the algebraic Bethe ansatz for the new models. 

In section 4 we turn to the supersymmetric case and add multiplicity to integrable models 
associated with the su(m| n) super algebras. The form of the algebraic Bethe ansatz for these 
models is similar to those of the usual su{m\n) chain [H], but with modifications due to the 
additional degrees of freedom. It will turn out that this alters the spectrum of the monodromy 
matrix and leads to a degeneracy of eigenstates just as in the purely bosonic case. 



2 A„^_i-type Integrable Models With Multiplicity 

As a warm-up, and to establish notation, we shall consider first bosonic models similar to those 
in [7], before we turn to the supersymmetric case in section HI Let {Ej : / = 0, . . . , m — 1} be the 
standard basis oi V = C"\ {fi^} the dual basis of V*, and E/ = Ej ® e V ® V* = End{V). 
Suppose R°{u) G End(y(S)V") is an i?-matrix obeying the usual Yang-Baxter equation with spectral 
parametei0 

Rl,{u - v)RUu - w)Rl^{v -w) = R°,{v - w)Rl^{u - w)Rl^{u - v), (2.2) 
which, moreover, is of the particular form 

m—1 

^°(^) = H {R'i'jW E/ ® E/ + (1 - 6jj)R'/j{X) E/ ® E/) . (2.3) 

I,J=0 

That is, in each non-zero element of the i?-matrix, the set of 'in' indices is the same as the set of 
'out' indices. The usual trigonometic i?-matrix of A^^i |Sj is of this type, and was the particular 
case considered in [7]. From the point of view of statistical mechanics this property originates in 
the ice rule obeyed by the Boltzmann weights. 

Given any such i?-matrix, the "multiplicity" models of [7] are defined by replacing each in- 
dividual basis state by a collection of states. Thus, for every Ej, introduce new basis states 
{caj '■ ai G ■£^i}i where ^ is an index set of cardinality rij G N. Let {u""' : a/ G s^i} be the dual 

^Wc follow the standard notation of [10], so that for example Rij,{u) denotes the element of End(y)'^'^ 

m—1 

^ R'j''l{u)Ei'®\®Ek''. (2.1) 

I,J,K,L=0 
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basis, write Caj'^-' = e^^ ® u^-', and then define 

1 

R{X) =J2Y1 (^Vj(A) eb/' ® e,/^ + (1 - (5,j)i?^/j(A) e,/^ ® e^/^) . (2.4) 

/,J=0 ajGii// 

One way to see that this new i?-matrix also obeys the Yang-Baxter equation is to note that R°{X) 
may be written purely in terms of the combination Ej (g) (no sum on /), with the tensor factors 
inserted in appropriate slots. The new i?-matrix is obtained by performing the replacement 

Ei^n^ J2 (iai ® i^"' , (2.5) 

SO it suffices to observe that these objects obey the same rule under internal contraction: on the 
one hand 

Ej {n^,Ej)(^n-^ = SijEj(^n^, (2.6) 

and on the other, since the index sets j^i are disjoint, 

J2 ea, U'^', 5^ \ ® a;^^ = 5^ <5a,6,e,, ® a;^'^ = e^, ® uj"" . (2.7) 

We consider a spin chain of length po, and take as our Lax operator (see e.g. [11] for a review) 

L^M = Rax{^-Xl), (2.8) 

where x G {1,2, . . . ,po} denotes the space of the xth particle on the chain, while the index a 
denotes the auxiliary space. The A° appearing in the arguments of the /^-matrices parametrize 
the inhomogeneities at each position of the chain. From the Lax operator, we construct the 
monodromy matrix 

{fi)---L^M. (2.9) 

The transfer matrix is the partial trace of the monodromy matrix: r(/i) = tr^, TQ,(/i). Notice that 
the real difference between the model with multiplicity and the usual model the lies here: the 
monodromy matrix makes explicit reference to the indices a/, and the notational trick used above 
to recast the new i?-matrix in the same form as the original one cannot be used here to simply 
reduce the new model to the old one. 

3 Algebraic Bethe Ansatz 

The algebraic Bethe ansatz for these models has the same structure as the one developed by 
Maassarani for his models in [7j. In fact, the models of [7] are subsumed by the more general 
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framework here, though we have not included the parameters Maassarani denotes as Xq,.^^ , which 
in our notation would be Xajbj- Adding these parameters is not difficult, but we will work with 
them set equal to unity. 

To proceed with the algebraic Bethe ansatz, choose a pseudo-vaccum: 

l^o) = |ao---ao), (3.1) 

where ao G ^ and the capital Latin indices now split into / = (0,2). Let a/(A) = i?{{(A), 
= RjjW and cij{X) = Rjj{\) for I ^ J. The choice of pseudo- vacuum singles out of 
various elements of the monodromy matrix. Let A{\) = T^^{X) and Baj{X) = r^"|'(A). It is not 
hard to check that on the pseudo-vaccum we have 

po 

A(/i)|fio) = llaoifx-Xlmo), 

x=l 
Po 

T::\no) = llko{^^-\l)\no). (3.2) 

x=l 

The only other operators that do not annihilate the state are 

Po 

Bao{\)\^o) = JJao(/i- A°.)|aoao- --ao), 

x=l 

PO / PO \ /x-1 \ X 

5a,(A)|f^o) = I n '^o(/W- A°) j Cio(/i- A°) I Jj6io(/i- Aj) j |ao---''^---ao). 

x=\ \z=x^\ J \y=l J 

(3.3) 

The second of the above two relations is what we will use to construct the other eigenstates of the 
transfer matrix. These states will be of the form 

m— 1 

1^1) = E E F'=^^^'=^-^'="i?c,,(A})Se,,(A^)---Sc,,^(AiJ|f]o), (3.4) 

where F is to be determined. We break the analysis of the action of the transfer matrix up into 
m — 1 levels and at each level aside from the final one, the procedure remains essentially the same. 
First, work out the action of T^°(/i) on where Oq 7^ ciq. These elements of the monodromy 

matrix annihilate |f2o), but it turns out that they preserve states of the form |aj^ ■ ■ ■ cLi^,^) when 
Pi = Po- This is an effect of the multiplicity that does not occur in the usual models. Explicitly 
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we have 

m— 1 

Po 

= WbQi{ii-\l)\a^^---ai^^). 

x=l 

As a result, the action of on is 

r,7(/x)|vl/i) =5„„, J]6o.(^-A°.)|*i). (3.5) 

X' = l 

The action of remaining components of monodromy matrix can be dealt with using the RTT- 
equations, 

R^p{\ - ^x)T^{\)Tp{^i) = T^{fi)T^{X)R^^{X - /i), (3.6) 
which, with the matrix indices displayed explicitly, read 

E E (A - a^t;;;^ m:::; (/^) = E E ^-t (/^) A (A)<Cr (a - /.). 

Mi,M2 eMiS^/A/i Mi,M2 eu^^-sfhl^ 

(3.7) 

By specifying the various free indices, we may write down the particular relations which interest 
us: 

A{^)B^.W = ' ^W w^if^) - " ^I bM'^W (3.8) 

bio{X - jJL) bio{X - /i) 

Ignoring contributions from the second term on the right-hand-side of either equation, we proceed 
by computing first 

^(/^)i*i) = n «o(/^ - A°) n i^i). (3.10) 

x=i y=i^o(A,y-^j 
where we are forced to make the simplification bko = &o in order for to be an eigenstate of 
j4(/i). Now compute the action of T^' on 

po Pl m— 1 



X s;,^ (A}) ■ ■ ■ Bj,^^ (xijino) (Ti)::;S;:::5^ (^)f-i--.i , (3.11) 



where the level 1 monodromy matrix is defined as 

T„\(^) = L,,p,(/i - AJJ ■ ■ • - A}), (3.12) 

and ai refers to the level 1 auxiliary space whose components only run over the indices in the 
sets . Clearly, we have now come around to an analogous problem one level up. The 

undetermined coefficients F appearing in are in fact, states of a level 1 spin chain acted on by 
the level 1 monodromy matrix. We wish to choose these to be eigenstates of the diagonal blocks 
of the level 1 monodromy matrix: 

(Ti)^:(^)F=(Ai)::(/i)F, (3.13) 

where (A^)^^(;u) is the level 1 eigenvalue. Given this, the expression for the action of r^*|\I^i) 
simplifies to 

PO pi ^ 

= n - A°) n . . __xn (^')SMI^i)- (3.14) 

x=i y=i 

So we have deduced the eigenvalue up to the next level in the analysis. It is 

PO PO PI / \ 1 _ \ 

A(/^) = ^popiK-i)n^o(/^-A°)+n«o(/i-A°)n ' 



x=l x=l 

PO Pi H 

+ nM.-A»)ni;^^^A-w. 



x=l 

where = e^/ (A^)2- '^^^ factor of rzo — 1 appearing in the first term arises from the part 

of the trace that sums over all the states in j^q except for oq. The higher levels have the same 
structure, so the higher level eigenvalues are 

Pk Pk Pfe+l ^ (X^+l _ n) 

a'(ao = w-.i(^^ - 1) n ^^-(^ - + n ^^'(^ - n _ 

x=l x=l y=l ^^^-^y 

Pk Pk+l ^ 

+ riMf-AarirT— i^A'«M^ (3^15) 

x=i y=i "^^^^ ^y I 

Here k runs over 0, . . . , m — 2. We have also defined 6fc(/i) = bki{fi) for alH = 1, . . . , m — 1. Also 
note that throughout we have left implicit the dependencies on the inhomogeneities {Aj^, . . . , A^'^j 
that appear at each level k = 0,...,m — 2. The final eigenvalue A"^~^ appears when k = m — 2. 
Examining the i?TT-equations for this final level: 
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where all the indices above take values in s^m-x- The level m — 2 creation operator is given by 
(S"^~^)c = where 0^-2 ^ is the chosen level m — 2 pseudo-vacuum state. It is 

not hard to see that for e, /, a, c G ^m-i we have 

KiW = am^i{X)5[5^, = a^_i(A)P:/, (3.16) 
where P is the permutation operator. So the equation simplifies to 

- ^3";|^:^| (i^"^-^)aM(T-^)^(A). (3.17) 

We can deduce that the level m — 1 monodromy matrix is built out of permutation operators: 

Pm-l 

{T--')M = P^,_,---P^^ n «(/^-^r')- (3.18) 

x=l 

The product of permutation operators yields the unit-shift operator for a chain of length Pm-i- 
The eigenvalues of this operator are roots of unity, and thus, the eigenvalues of the level m — 1 
transfer matrix are of the form of a root of unity times the product Hx^i^ '^(a* ~ A^~^). 

The Bethe equations for this chain are derived by requiring that the residues around the various 
A-parameters vanish. Thus for A; = 0, . . . , m — 2 we have 

^ fx «.(A^^ - A^) «.(Ar^ - A^^)6.(A^^ - A,^^^) 
M&.(A^^-A^) Jl-L f>,(A^^+i-A^i) 

where we have assumed that 6fc(/i) — when ~ 0. Substituting in the expression for A'^+^(A^"*"^) 
and rearranging yields 

11 h. . J\k+2 _ \k+l\ LL n,(\k+l - \k+l\h.(\k+l - \k+l\ 11 



\ &fc+i(A,^+2 - A^i) J-jL «^(A^+i _ A^i)6,(A^i - A^-+i) 11 a,{X^+i - A^) a,+i(0) ' 
y¥=z 

(3.19) 



4 A^_i\n-i Models 

We now turn to our main topic of interest, namely the introduction of multiplicity to supersym- 
metric chains. Suppose henceforth that y is a Z2-graded vector space and that our basis is chosen 
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such that each basis vector Ej has a definite grade |/| G {0, 1}. Because vectors of grade (respec- 
tively 1) are taken to obey bose (fermi) exchange statistics, the tensor product becomes braided 
(see for example [121 [131 E]) in a mild fashion: 

(1 ® Ej){Ej ® 1) = (-l)l^ll'^l^j ® Ej. (4.1) 

This braiding means that there is some subtlety when taking tensor products of elements of 
End(y), essentially because the isomorphism between End(y) ® End(y) and End(y ® V) is no 
longer quite trivial. Let us regard Ej'^ = Ej ® ^l'-' ^ End{V) as acting on V* from the rig htS 

E/ : n^' — ^ n^^E/ = 6i^n^. (4.2) 

We then define ®s by the demand that, acting from the right on V* ^V*, E/"^ ®s Ej^ send 

Noting that Ej^ is of grade |/| + \K\ mod 2, we have 

(fi^^®fi^)(^,^®£;/) = (-l)l^l(l™)(n^E,^®fi^E/) (4.4) 

= (-l)\Jm+W)6f6^ Q^'^Q'^. (4.5) 

and therefore 

Ej^ ®, E/ := (-1)1-^1(1^1+1^1)^^^ ® Ej^. (4.6) 

The definition extends naturally to more copies of End(V^), by including signs as needed to ensure 
that the action from the right on V* ^ ■ ■ ■ V* is correct. ®s is often called the "graded" 
or "supersymmetric" tensor product, although with the present conventions this is somewhat 
misleading because the grading is already built into (8>, as in fl4.ll) . 

Given these preliminaries, it is possible to proceed much as in the bosonic case, being careful to 
allow for the extra signs in (14. 6p . Consider an i?-matrix obeying the graded Yang-Baxter equatioijf 

Rl,{u - v)R'',,{u - w)Rl^{v -w) = R°,{v - w)R1,{u - w)Rl^{u - v), (4.7) 

which is of the particular form 

m+n—l 

R"W= E {rijWE/^sE/ + {-iy'^^'^{l-6ij)tij{X)EjJ(^sE/). (4.8) 

LJ=0 



^This will allow our conventions to match those of [8|; we could instead consider the action from the left on V 
but the resulting definition (|4.6p would be rather different. 

•^The advantage of using (E)s appears here: this is an equation in End{V (E}V (E}V) = End(y)(8)sEnd(y)(8)sEnd(y), 

and in component form the multiplication is straightforward matrix multiplication. The extra signs are wrapped 
up in the definition of e.g. R12 = R^s 1. 
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One class of examples are the su{m\n) i?- matrices of [131 E], whose Boltzmann weights are 

rn{X) = aj{X) = ^- , _ ' (4.9) 

a(l - e'^^) 

ruW = b{X) = I^J (4.10) 

tu{X) = c{X) = , />J (4.11) 

9 e 



t,,(A) = rf(A) = J<J. (4.12) 

The algebraic Bethe ansatz for models built from this i?-matrix (including integrable impurities 
involving dual representations) was performed in [8]. 

Now we introduce multiplicity to the models. Once more, introduce the index sets {-2^}^""^, 
and the vector and dual vector bases {ea/}^g^7=o' {'^'^^}a^e^i"=o- These have precisely the same 
properties as before, and again the modified i?-matrix 

m— 1 

^(^) = E Z (^^■^(^) ®^ + - ^/J)^/J(A) ei>r ®s . (4.13) 

/,J=0 a/Gii// 

satisfies the Yang-Baxter equation. We choose the Lax operators 

L^M = R^^{^-\1) (4.14) 
and the monodromy matrix for an inhomogeneous chain of length po is then 

TaifJ-) = -^«po (/^)-^a(po-l) {f^)---LM- (4.15) 

The transfer matrix is obtained by taking the supertrace in the auxiliary space a: 

r(/i) = str„r.(/i) = J2J2 (-)'"^a7(/^)- (4.16) 

It follows from the Yang-Baxter equation that T(^) satisfies the RTT relations 

Ri2{\ - fi)T,{X)T2{fi) = T2{^)Ti{X)Ru{X - fx), (4.17) 

With all indices displayed explicitly, these read (some care is needed with the sign in T2(/u)Ti(A) = 
(1 ®s T{fi)){T{X) ®s 1) on the right-hand side here) 

(_l)|L|(|J|+|M|)^/^i.^(^ _ ^)T^j^^^^X)TfL^^^^) = (-l)\-mi\MP\)T^j^^^^)T^i^^^X)R\jQ^{X - /.). 

(4.18) 
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The algebraic Bethe ansatz for the supersymmetric XXZ models is very similar to that of the 
ordinary XXZ models, the main difference being a matter of including the signs in the equation 
above. This means that the algebraic Bethe ansatz for Am\n models will be similar to that of the 
previous purely bosonic models. Just as before we choose a pseudo-vaccuum: 

\no) = \ao---ao), (4.19) 

where ao E s^o and the capital Latin indices now split into / = {0,i). The elements of the 
monodromy matrix that this singles out are A{X) = T^°{X) and i?5j(A) = TI^°{X). We are interested 
in the action of the following operators on the pseudo-vacuum 

PO 

A{f,)\no) = Haoifi-Xlmo), 

x=l 

PO 



t::W = llb{^l-x'^)\no), 



x=l 

PO / PO \ /x-1 



x=l \z=x+l / \y=l 



(4.20) 



The other eigenstates of the transfer matrix will be linear combinations of the states 



m—l 



l^i) = E E ^^'=^^'=^-^'^-^c,,(A})i?e,,(A^)---5e,,^(A;j|fio), (4.21) 

i=i cfe^ G£4i v.cfcpj e.64p-^ 

where F will turn out to be a chain at the next level of nesting. Recall from before that the 
multiplicity leads to contributions that do not appear in the standard models. Namely, 

pi 

rC(/x)|^i) = 6,,,, n - (4.22) 

x=l 

We may read off the relations that we need from the RTT equations 

T^if^B,.w = E E (-)i*i(i-H.i)H-ioi(i.iH-bi) ^a:y^-^ 

r,s=l er€^r 

- (-)'°ll^l^l^'l(l°l^l^l^||^5.(/i)r,^(A). (4.24) 

Observe that the i?- matrix appearing in the second equation is the one that describes the Am-i+\o\\n- 
model. Since each additional level (until the final one) is identical in structure to the initial level. 



10 



we deduce that the general eigenvalue at each level is 

Pk Pk Pk+1 n (\k+l _ 

x=l x=l y=l '^'^^y 

Pk Pk+l 

+ nK^-^^)n ^(. A^H-i) ^'"^M> (4.25) 

x=i y=i ^ y ' 

for /c = 0,...,Tn + n — 2. The level A; + 1 transfer matrix is defined as 

= ^ ^ (_)E^iT'ioiir|,+E^=i{k.iir.i+Nis,i) 

i,ri,...,rpj._i=A:+l aiG.c/i 

(4.26) 

We can derive this transfer matrix from a suitably defined level k + 1 monodromy matrix 

ri+\(/x) = L^,,,pUf^ - A5+\) ■ ■ ■ L^,,Af^ - A^^). (4.27) 

The products in the level k + 1 auxiliary space whose components only run over the indices in the 
sets {-i^^l^fc+T^ supertrace are defined to account for the signs appearing in fl4.26p . 

The F coefficients at level k are chosen to diagonalize the level A; + 1 transfer matrix 

T''+\fi)F'' = A'+\ij)F\ (4.28) 

where A^"'"^(yu) is the level k + 1 eigenvalue. 



At the final level the i?TT-equations are the same as those appearing in the bosonic case 

b^_i{lj, - A) 



(T"-^)a(/^)(i?""')c(A) = ^/^ (i?-^);(A)(r-^)^(/i) 



Cm-l{fJ' A) , , ^_2sft. 



bjn-l{^ - A) 

The indices take values in s^m-i and the level m — 2 creation operator is (i?™~^)c = (T™"^)"'"-^, 
where am-2 ^ is the chosen level m — 2 pseudo-vacuum state. Just as in the bosonic case, 

the level m — 1 monodromy matrix will turn out to be proportional to a product of permutation 
operators yielding the unit shift operator. Thus A™~^(^) will be a root of unity. 
For k = 0, . . . ,n + m — 2 the Bethe equations are 

a.-+i(Ar - A^^) ^r^ a.+i(Ar^ - A,^+^)b(A^^ - A^^) -p^ 6(Ar^ - A^) _ 
11 6(A^+2_^fc+i) 11 a,(A^^+i-A^i)6(A^i-A^'+i) H «,(A,^+i - A^) ' 

(4.29) 
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where we have used the exphcit formulae for the functions ai{^) and b{fi) to cancel off some 
signs and simplify the expression a little. Notice that these are essentially the same as the Bethe 
equations for the bosonic case with a/(0) = 1. 

The spin chain Hamiltonian for these models is the logarithmic derivative of the transfer matrix 
with respect to the spectral parameter A. The calculation is simpler when one assumes that all 
the level-0 inhomogeneities vanish since the Hamiltonian is then calculable as the derivative of the 
i?-matrix. Using the explicit i?- matrix given by the coefficient functions (14. 9 p - ( I4.12p . we have 

dR{X) 



H = V- 



dX 



/ ^ — ' ^ — ' \ i — o , 



_ W\J\e bj ^ a; 



(tS) E E ®^ -^/^ + (r^) E E e./^ (4.30) 



I<J ai^s^l ^ ' I>J aje.e/j 



where V is the graded permutation operator 

m+n~l 

r= E (-)"""^'^/'®^^^/'- (4-31) 

I,J=1 ajG.a'i 
bj£.e/j 

The results of the algebraic Bethe ansatz allow us to compute the energy spectrum corresponding 
to this Hamiltonian. If we make the drastic simplification that po > pk for all A; > 0, then we can 
arrive at a closed form expression. By taking the logarithmic derivative of A°(A) we find that if 
the level-1 states are bosonic |A; = 1| =0, then 

^ sinh 7 / . X 

E=} ^- ( -, (4.32) 

^smh (A^ + 7) smh A^ 

where we have let q = e~'^. The answer differs a bit if the level-1 states are fermionic |A; = 1| = 1: 

cosh 7 ^coshA^ 
E = po-p,- Po—r^ + E — rrf ■ 4.33 
smh 7 ^-^ smh A,, 

y=i y 

5 Conclusions 

In this paper we extended the multiplicity Am models of [7] to the supersymmetric case. We 
used the nested algebraic Bethe ansatz to find the eigenvalues of the transfer matrix and the 
corresponding Bethe equations. Using these results, we computed the Hamiltonian and its energy 
levels for the specific multiplicity model associated with the supersymmetric XXZ model. 
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While the multiphcity of the model can be hidden notationally in the description of the model's 
i?-matrix, the addition of multiplicity is reflected in the form of the nested ABA equations. The 
final level of the nesting becomes more non-trivial, as in fl3.18p . and there are also additional terms 
proportional to 5^^^^+! at every other level. The multiplicity also contributes to the energies in 
the corresponding spin chain when po = pi or (some of) the level-0 inhomogeneities are non-zero. 

Since the XX model is a building block for more complicated theories, notably the Hubbard 
model, it would be interesting to investigate theories constructed by coupling together supersym- 
metric multiplicity models. Generalizing the Hubbard model while maintaining integrability has 
proven difficult [13, [THl [13 [TSl HSl 120] • We hope that the work presented here will shed some light 
on this problem. 

Acknowledgements 

We are grateful to Jonathan Evans for helpful conversations. The research of C.A.S.Y. was 
funded by the Lever hulme trust. 

References 

[1] Z. Maassarani and P. Mathieu, The su{n) XX Model, Nucl. Phys. B 517, 395 (1998) 
[arXiv:cond-mat /9709163] . 

[2] Z. Maassarani, The SU{n) Hubbard Model, Phys. Lett. A 239 (1998) 187-190. 

[3] R. Yue, R. Sasaki Lax pair for SU{n) Hubbard Model [arXiv:cond-mat/9801193v2] 

[4] Z. Maassarani Exact integrability of the su{n) Hubbard model Mod. Phys. Lett. B 12 51 
(1998) 

[5] M. J. Martins, On the integrability of the SU{N) Hubbard model, Phys. Lett. A 247 (1998) 
218 [arXiv:cond-mat/9710049vl] 

[6] Z. Maassarani, The XXC models Phys. Lett. A 244 (1998) 160 [arXiv:solv-int/9712008] 

[7] Z. Maassarani, Multiplicity Am Models, Eur. Phys. J. B vol. 7, 627 (1999) 
[arXiv:solv-int /9805009] . 

[8] G. A. P. Ribeiro and M. J. Martins, Algebraic Bethe ansatz for an integrable Uq[Sl{n\my\ 
vertex model with mixed representations, Nucl. Phys. B 738, 391 (2006) 
[arXiv:nlin.si/0512035]. 



13 



[9] M. Jimbo Quantum R matrix for the generalized Toda system Comm. Math. Phys. 102 
(1986) 537-547 

[10] L. D. Faddeev, N. Y. Reshetikhin and L. A. Takhtajan, Quantization of Lie Groups and Lie 
Algebras, Leningrad Math. J. 1, 193 (1990) [Alg. Anal. 1, 178 (1989)]. 

[11] A. Kundu, Quantum Integrable Systems: Construction, Solution, Algebraic Aspect, 
[arXiv:hep-th/9612046]. 

[12] S. Majid, Foundations of quantum group theory, Cambridge Univ. Press (1995) 

[13] Y. Z. Zhang, Comments on the Drinfeld realization of the quantum affine superalgebra 
Uq{gl{m\nY^^) and its Hopf algebra structure, J. Phys. A: Mathematical and General 30, 
8325 (1997) 

[arXiv:math.qa/9809156]. 

[14] Y. Z. Zhang and M. D. Gould, Quasi- H op f Superalgebras and Elliptic Quantum Supergroups, 
J. Math. Phys. 40, 5264 (1999) 
[arXiv:math.qa/9809156]. 

[15] P. Coleman, C. Pepin, J. Hopkinson, Supersymmetric Hubbard Operators, 2001, Phys. Rev. 
B 63, 140411 

[arXiv:cond-mat/0007495]. P. Coleman, J. Hopkinson, Atomic Model of Susy Hubbard Oper- 
ators, 2003, Phys. Rev. B 67, 085110 
[arXiv:cond-mat /0202060] . 

[16] N. Beisert, The su{2\2) dynamic S-matrix, 
[arXiv:hep-th/0511082]. 

[17] F. H. L. Essler, V. E. Korepin, K. Schoutens, New Exactly Solvable Model of Strongly Corre- 
lated Electrons Motivated by High-Tc Superconductivity, Phys. Rev. Lett. 68 (1992) 2960. 

[18] A. J. Bracken, M. D. Gould, J. R. Links, Y-Z. Zhang, A New Supersymmetric and Exactly 
Solvable Model of Correlated Electrons, Phys. Rev. Lett. 74 (1995) 2768 
[arxiv:cond-mat / 9410026] . 

[19] J. M. Drummond, G. Feverati, L. Frappat and E. Ragoucy, Super-Hubbard models and appli- 
cations, JHEP 0705, 008 (2007) [arXiv:hep-th/0703078]. 

[20] M. J. Martins and P. B. Ramos, Solution of a Supersymmetric Model of Correlated Electrons, 
Phys. Rev. B 56 (1997) 6376 [arXiv:hep-th/9704152]. 



14 



